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Correlation and time delays of the X-ray and optical 
emission of the Seyfert Galaxy NGC 3783 
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ABSTRACT 

We present simultaneous X-ray and optical B and V band light curves of the Seyfert 
Galaxy NGC 3783 spanning 2 years. The flux in all bands is highly variable and the 
fluctuations are signiflcantly correlated. As shown before by Stirpe et al. the optical 
bands vary simultaneously, with a delay of less than 1.5 days but both B and V 
bands lag the X-ray fluctuations by 3-9 days. This delay points at optical variability 
produced by X-ray reprocessing and the value of the lag places the reprocessor close 
to the broad line region. A power spectrum analysis of the light curve, however, shows 
that the X-ray variability has a power law shape bending to a steeper slope at a time- 
scale ~ 2.9 days while the variability amplitude in the optical bands continues to grow 
towards the longest time-scale covered, ~ 300 days. We show that the power spectra 
together with the small value of the time delay is inconsistent with a picture where all 
the optical variability is produced by X-ray reprocessing, though the small amplitude, 
rapid optical fluctuations might be produced in this way. We detect larger variability 
amplitudes on long time-scales in the optical bands than in the X-rays. This behaviour 
adds to similar results recently obtained for at least three other AGN and indicates 
a separate source of long term optical variability, possibly accretion rate or thermal 
fluctuations in the optically emitting accretion disc. 

Key Vifords: Galaxies: active 



1 INTRODUCTION 

The optical continuum emission in Active Galactic Nuclei 
(AGN) almost certainly originates as thermal emission from 
the accretion disc (Koratkar & Blacs 1999) but the origin of 
the X-ray emission in these objects is less clear. The accre- 
tion disc is expected to be too cool to produce X-rays and 
the spectrum in this band is non-thermal, so the existence of 
a 'corona' has been suggested to produce X-ray emission by 
Compton up-scattering optical/UV photons to higher en- 
ergies. As it has not been established how this corona is 
formed, its location and relation to the accretion disc are un- 
certain, although variability and energetics arguments place 
it in the innermost regions of the AGN. One way to probe the 
connection between the accretion disc and corona is to track 
the variability in X-ray and optical bands simultaneously, to 
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determine their degree of correlation and the relative loca- 
tion of their emitting regions by measuring delays between 
the fluctuations in both bands. 

If the X-rays are produced by a corona in the vicinity of 
an optically thick accretion disc, it is expected that at least 
some of the X-ray flux will be intercepted by the disc and be 
reprocessed thermally into optical emission. The short vari- 
ability time-scales observed in optical bands and the short 
delays with respect to the X-ray variations suggest that the 
optical variabili ty arises from repro cessing of the highly vari- 
able X-ray flux ( Krolik et al.lll99lh . The radial temperature 
profile for a standard accretion disc, T{R) oc R~^^^, places 
the main optical emitting regions at different radii in the 
disc, leading to longer light travel times for longer wave- 
length optical emission. The wavelength dependent delays 
predicted by this model are in agreement with the obser- 
vations of several AGN (IWanders et al ] |l997l : ICollier" et al.l 
l200ll : ICackett et akllioOTl V 
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There is an increasing number of objects, however, 
whose X-ray/optical behaviour is in conflict with a model 
where reprocessing is responsible for all the optical variabil- 
ity. Long time-scale optical fluctuations have higher ampli- 
tude than their X-ray counte rparts in at lea st three cases. 



NGC 5548 ( Uttlev et al 
[2OO8) and Mkn 



79 



20031'). MR 2 2 51-17 8 (jArevalo et al I 
Breedt et ahl 120091 '). X-rays have 
been seen to lag th e optical band s in e.g. NGC 4051 
f Shemm er et al ]\200± and Mkn 509 (jMarshall et al.ll2008l ) 
and Kaz anas fc NavakshinI l|200ll ') have modeled in detail the 
reprocessing in the case of NGC 3516 and shown that this 
process is unable to explain its optical variability. Repro- 
cessing is also challenged by energetics arguments in some 
AGN, where the luminosity in the optical bands largely ex- 
ceeds the luminosity in X-rays, making it impossible for the 
X-rays alone to pr oduce the obs erved large amplitude opti- 
cal flux variations (lGaskell|[2007l ). 

In I Arevalo et al.l (|2008l ) we proposed a composite model 
to explain the X-ray optical correlation in MR 2251-178. 
We allowed a small fraction of the optical emission to arise 
from reprocessing, in order to produce the observed small 
amplitude, rapid optical fluctuations, with a negligible de- 
lay with respect to the X-rays. At the same time, the long 
term variability in both optical and X-ray bands would be 
modulated by accretion rate fluctuations, allowing differ- 
ent amplitudes of variability in both bands on long and 
short time-scales. Predictions of this two-component sce- 
nario should be tested on AGN of different masses and ac- 
cretion rates, rh, since the location of the optical emitting 
region on a standard accretio n disc depends on these param- 
eters (e.g [Treves et al]|l988l ). Low m and high mass black 
holes should have a relatively cool disc, placing the optical 
emitting region at small radii, probably close to the X-ray 
emitting corona. Therefore, in these cases we should expect 
a stronger X-ray/optical correlation and variability on simi- 
lar time-scales. The opposite would happen for low mass and 
high accretion rate AGN. We are currently building a sam- 
ple of AGN monitored simultaneously in X-rays and optical 
bands, which covers a large range in mass and accretion rate 
to establish the connection between their optical and X-ray 
emitting regions. Here we present results on our third tar- 
get, NGC 3783. whic h has a black hole mass of 3 x lO^'M© 
l|Peterson et~aLll2004h a nd average accretio n rate m — 7% 
of the Eddington value (|Woo fc lJrrvll2002l ). 

The structure of the AGN in NGC 3783 has been stud- 
ied through several monitoring campaigns. Glasi (|1992| ) 
monitored NGC 3783 in the optical and infrared bands flnd- 
ing an ~ 80 day lag between 'V and K bands from a 15 year 
long light curve. This lag is interpreted as the light travel 
time between the central engine and the hot dust which 
reprocesses the optical/UV emission and is conflrmed and 
studied in detail by Lira et al. (in prep.). The measured lag 
is at precisely the value expected for the light travel time 
to the dust at the sublimation radius given the luminosity 
of NGC 3783, (see ISuganuma et al. 2006). The sparse sam- 
pling of the long light curves of [Glass (1992,), however, does 
not allow a measurement of lags between optical bands and 
the lack of X-ray coverage precludes a study of the origin 
of the optical variability. The structure of the broad line re- 
gion (BLR) has been studied through reverberation mapping 
campa igns lasting over 7 mont hs, reported in[Rei chert et al.. 
l| 19941 ) and lStirpe et all (|l994l '). These authors found almost 



simultaneous variability in the UV and optical continua with 
amplitude increasing with decreasing wavelength. They also 
find an ~ 8 day delay between the continuum and broad- 
line emission variations, establishing the distance between 
the BLR and the central continuum source. The optical/X- 
ray correlatio n in NGC 3783 on short timescales was in- 
vestigated by ISmith fc VaughanI ([2007,) from simultaneous 
X-ray and UV observations obtained with XMM-JVewton. 
The ~ 5 days long light curves in both bands show signif- 
icant variability but the data are not sufficient to detect a 
statistically significant correlation or lag. In this paper we 
study the optical/X-ray correlation on time-scales of days 
to hundreds of days using 2-year long, regularly sampled 
light curves with the aim of establishing the relation be- 
tween optical and X-ray emitting regions. The monitoring 
campaigns and data reduction are described in Sec. [2] We 
compare the variability properties of optical and X-ray light 
curves through the power spectrum in Sec. [3] and calculate 
their cross correlation and time lags in Sec. [l] We inter- 
pret the variability properties and lag time-scales in terms 
of physical models in Sec. [5] and summarise our conclusions 
in Sec. |6] 



2 DATA 

We have monitored NGC 3783 in the X-ray band using the 
Rossi X-ray Timing Explorer RXTE and in the B and Y 
optical bands with the 1.3m SMARTS telescope in Chile. 
The resulting light curves are shown in Fig. [1] Below we 
give a brief description of each observational campaign. 



2.1 X-ray monitoring with RXTE 

NGC 3783 was monitored by RXTE by taking approxi- 
mately 1 ks exposure snapshots at regular intervals. For the 
purposes of the present paper, we include X-ray data from 
2006 November 1 to 2008 August 31. The sampling pat- 
tern varied during this period: exposures were taken every 
two days from the beginning of the campaign to 2007 June 
29 then every four days until 2008 February 8. This was 
followed by an intensive sampling campaign when observa- 
tions where made three times per day from 2008 February 
8 to 2008 June 3, after this date the sampling returned to 
one observat ion every four days. The data were reduced as 
described in lArevalo et al.1 (|2008h and the light curve con- 
struction followed the procedure detailed in the same paper. 
In particular, the 2-10 keV flux was estimated by fitting 
an absorbed power law model to each snapshot spectrum, 
where the absorption column was fixed at the Galactic value, 
Nh = 8.7 X lO^^cm-^ (Lockman & Savage 1995). The long- 
term light curve, binned in 6 day bins is shown in the top 
panel in Fig. [1] 



2.2 B and V band monitoring with SMARTS 

The optical monitoring was performed using the ANDICAM 
instrument mounted on the 1.3m SMARTS telescope in 
Cerro Tololo, Chile. Two B and V band exposures of 
NGC 3783 were taken for each filter. Observations were 
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Figure 1. NGC 3783 light curves. From top to bottom: 2-10 keV X-rays, SMARTS data in the B band and V band. For clarity, the 
X-ray light curve has been binned in 6 day bins in this plot. The B and V fluxes were measured using psf fitting to the central luminosity 
peak of the galaxy images so they include some starlight contamination from the underlying galaxy. The intensive sampling in the optical 
light curves occurs at MJD-50000=4505-4580. 



made every four days between 2006 December 1 and 2008 
August 31 including a daily sampling period between 2008 
February 15 and April 28. The period of daily sampling in 
the optical was contained within the period of intensive X- 
ray monitoring. 

We performed psf photometry on NGC 3783 and four 
non-variable stars in the field of view using the IRAF task 
psf. We constructed relative flux light curves by dividing the 
flux measured for NGC 3783 nucleus by the sum of the refer- 
ence star fluxes. The errors were calculated by propagating 
the magnitude error produced by the psf task. We confirmed 
that the relative fiux of the reference stars was constant 
throughout the campaign. A B band image of NGC 3783 
and the reference stars is shown in Fig. [2] We note that the 
underlying galaxy fiux within the nuclear PSF is not sub- 
tracted from the AGN flux, so a relatively small amount of 
starlight contamination remains. 

We performed aperture photometry on the reference 
stars and used the zero-point correction factors calculated 
for the SMARTS telescope every photometric night, to- 
gether with the airmass of our observations, to calibrate 
their magnitudes. We converted these magnitudes into fiux 
densities at 5500 A for the V band and 4400A for B, as- 



suming a fiat spectrum {a — 0) within each band and used 
the total flux of the four stars to calibrate the relative-flux 
light curve of NGC 3783. The flnal light curve fluxes were 
corrected for foreground Galactic ex tinction using Ax = 
0.514 for B a n d Ax = 0.395 for V l|Schlegel et al.lll998l ). 
IWinkler et all (|l992l ) used several methods to estimate the 
extinction towards NGC 3783 and concluded that it can 
be largely or completely accounted for by the foreground 
Galactic value, we therefore made no further correction to 
the fluxes. 

The calibrated B and V band light curves are shown in 
the middle and bottom panels in Fig.[TJ respectively. The psf 
photometry method was used to include all the point source 
emission from the nucleus of the AGN while minimising the 
contribution from the host galaxy but it is likely that some 
contamination remains. We note however that in the B band 
light curve the flux varies from under 9 mjy to over 24 mjy, 
so the majority of the flux measured must be nuclear (as- 
suming that the star light is constant). In the V band the 
maximum to minimum variation is slightly smaller but still 
more than a factor of 2. This difference in variability ampli- 
tude can be caused by stronger star light contamination in 
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Figure 2. B band ANDICAM image taken with the 1.3m 
SMARTS telescope. NGC 3783 and the stars used as reference 
(SI to S4) are marked and labelled. The reference stars are not 
saturated and appear in all the campaign images. The N and E 
arrows are 1' long. 
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Figure 3. Intensive sampling NGC 3783 light curves. From top 
to bottom: 2-10 keV X-rays, SMARTS data in the B band and V 
band. The top panel shows the 3-times daily sampled light curve 
in crosses and these data binned in 1 day bins in circles. 



the V band but it is also possible that the AGN V emission 
is intrinsically less variable. 



3 POWER SPECTRUM 



The power density spectrum (PDS) can be used to 
quantify the variability amplitude as a function of the 
time-scale of the variations, or correspondingly, of their 
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Figure 4. Power spectra of NGC 3783 in the X-ray band (solid 
lines) and the optical B band (dotted lines). The dashed lines 
represent the best-fitting bending power model for each band. 
The apparent QPO in the X-ray light curve at / 5 X 10~^ Hz 
is not confirmed by the new intensively sampled data. The X- 
ray band shows clearly more variability power at high frequencies 
than the B band but at low frequencies (~ 5 X 10~* Hz) their 
power spectra intersect. 



Fourier frequency. The PDS is constructed through the 
modulus squared of the discrete Fourier Transform (DFT) 
iPress et al.ll 19921 ). For the normalisation used in our calcu- 
lations, the integral of the PDS over frequency equals the 
normalised variance of the Ught curve. 

For most AGN X-ray light curves, the PDS has a power 
law shape of slope ~ — 1 bending to a steeper slope at high 
frequencies (e.g. Summons et al. in prep., McHar dv et al.l 
12004 l2005l ). which is similar to the PDS fou nd in stellar 
mass black hole binaries in the soft state (see IUttlevll2007l . 
for a review). It is customary to multiply the variability 
power by frequency when plotting the PDS, to highlight 
deviations of the power law slope from -1 as, in this case, 
the low frequency part of the PDS appears approximately 
flat and the breaks are more noticeable. We use this standard 
of pres entation in Fig. |4] 

In ISummons et air feOO?') we show the X-ray PDS of 
NGC 3783 and explore the significance of an apparent quasi- 
periodic oscillation (QPO) at a frequency of ~ 5 x 10~''Hz. 
The new, intensively sampled light curve obtained for this 
object, shown in Fig. |31 with a sampling rate of three times 
daily and a length of four months, covers the frequency 
range 10"*^ — 3.4 x 10~^ Hz. This range covers the time- 
scales corresponding to the peak frequency of the possible 
QPO very well and allows us to test its significance conclu- 
sively. We calculated the PDS using the long term RXTE 
light cur ves and short term XM M-Newton light curves dis- 
cussed in [Summons et all l|2007l ) and added the new inten- 
sive RXTE data. The resulting PDS is shown in solid lines 
in Fig. 31 where the segments correspond to the different 
X-ray light curves used. 

We fitted a bending power law model defined as 



Pif) 



Af- 



1 + if/fb)""' 



(1) 



to the PDS using the Monte Carlo fitting technique psresp 
of lUttlev et al.l l|2002l ) . The low- frequency slope aL , the high 
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Figure 5. Confidence contours for the PDS parameters an and 
/;,. The sohd and dashed lines represent the 67% and 90% proba- 
bility contours, respectively. The area on the left of the open solid 
line correspond to the best-fitting region for the B band data and 
the closed curve on the right contains the best-fitting parameter 
values for the X-ray data. The parameter spaces allowed for each 
of the bands do not overlap showing that the PDSs are signifi- 
cantly different. 

frequency slope an, the bend frequency and the normal- 
isation A were allowed to vary. The best-fitting parameters 
are ql = 0.8, an = 2.2, = 5.8 x 10~^ Hz, co nsistent 
with the values obtained bv ISummons et ahl (|2007l ) for the 
same bending power law model. The corresponding model 
is shown by the dashed line in Fig. The simple bending 
power law provided an excellent fit to the new data (90% 
acceptance probability), making the possible QPO feature 
unnecessary. 

We also computed the PDS of the B band data, shown 
by the dotted line in Fig. |4l The long term light curve was 
used to constrain the PDS at frequencies 3 x 10~* — 2 x 10~^ 
Hz and the intensive light curve covered the range 2 x 10~ — 
7 X 10"'' Hz. We fixed ql = 0.8, i.e. the best-fitting value 
found for the X-ray PDS, for direct comparison with those 
data. Fig. [5] shows the 66% (solid lines) and 90% (dashed 
lines) confidence contours for an and of the optical and 
X-ray data. The parameter space regions do not overlap at 
90% significance, showing that the power spectra are signif- 
icantly different (i.e. they cannot have the same an and 
values simultaneously). The best-fitting bend frequency is 
below the lowest frequency probed and the B band PDS is 
consistent with a single steep power law model. 

4 X-RAY/OPTICAL CROSS CORRELATION 

We calculated the cross correlation between the X-ray and 
optical light curves using the discrete c orrela tion function 
(DCF) described in lEdelson & Krolik' (l988') and the z- 
transformed DCF of [Alexander, (1997). These cross corre- 
lation functions measure the degree of correlation between 
the two light curves as a function of r, the displacement of 
one of the light curves on the time axis. A positive value of 
the T will indicate that X-rays lead the optical fiuctuations. 

We first computed the long term DCF between X-ray 
and B and V light curves. The intensive sections were resam- 



pled by taking 1 point every 4 days in the optical bands and 
1 point every two days in the X-rays, to match the sampling 
of the rest of the corresponding light curve. As the rest of 
the X-ray light curve was sampled partly every 2 and partly 
every 4 days, we binned the re-sampled data into 4 day bins 
to obtain uniform sampling and probe similar time-scales 
throughout the light curve. We then computed the DCF for 
each year-long segment and combined the resulting func- 
tions, weighting each segment by the number of points in 
each lag bin. The B and V band light curves are very simi- 
lar so the DCF between either of them and the X-rays are 
virtually identical. The long term B vs X-ray DCF is shown 
in Fig.|6] As a separate check we also used the complete light 
curves (both segments together) resampled and binned into 
4-day bins as above and computed the z-DCF. The results 
were consistent with the DCF results. 

We estimated the significance of the DCF by calcu- 
lating the DCF between the observed optical light curve 
and 2000 simulated, uncorrelated X-ray light curves. The 
simulated X-ray light curves were generated following the 
PDS form and parameters describ ed in Sec. [3] using the 
method of lTimmer fc Konid l|l995h . and resampled and re- 
binned exactly as the real data. This method used to pro- 
duce the simulated light curves uses random phases for the 
power spectrum components, while in the real data there 
might be coupling between these phases, that produces for 
example the rms-flux re lation in the X-ray light curves 
(|Uttlev fc McHardvl[200ll '). We therefore note that we did 
not take into account the non-linearity of the light curves 
when estimating the significance of the DCF, which might 
have an efi'ect on the significance levels. For each simulation, 
the DCF was calculated for each of the two segments and 
then combined, just as was done with the real light curves. 
The median, 95% and 99% limits of the distribution of simu- 
lated DCF values are plotted in Fig.[6]in dot-dashed, dotted 
and solid lines, respectively, showing that the central peak 
measured between the real X-ray and B band light curve 
is higher than more than 99% of simulated data. No other 
peaks or dips in the DCF between the X-ray and either B 
or V band appeared significant within r ~ 200 days. We 
repeated this procedure applying the z-DCF to the com- 
plete real and simulated light curves in a consistent way. 
The central peak in the z-DCF between X-ray and B bands 
was significant at the 98.4% level and between X-ray and V 
bands at the 97.5% level. 

The short time-scale correlation was probed in more 
detail using the intensive-monitoring light curves shown in 
Fig. These include daily sampled optical data for 2.5 
months and X-rays sampled three times per day for four 
months. As shown above in Sec.|3l the X-ray light curve has 
considerably more variability power than the optical bands 
on these shorter time-scales. Fortunately, as the sampling of 
the X-ray band is more intensive we are able to cancel part 
of the high frequency variability in this band by binning the 
light curve in 1 day bins. In this way, the correlation between 
fiuctuations on time-scales longer than a day in both bands 
becomes clearer. We calculated the DCF between the 1-day 
binned X-rays and the V and B band intensive light curves. 
The significance of the correlation was estimated as above 
and the central peak of the DCF was higher than 99% of the 
2000 simulated DCFs for all light curve pairs. The intensive- 
sampling DCF is shown in Fig. [7] together with the signifi- 
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Figure 6. The DCF derived using all of the available long 
timescale X-ray and B-band data shown in Fig. [T] is plotted as 
black crosses, positive lag values correspond to X-rays leading. 
The solid lines represent the 99% extremes of the distribution 
of simulated X-ray light curves when correlated with the real B 
band light curve, the dashed lines represent the 95% extremes 
of this distribution. The correlation peak around time lag=0 is 
significant above the 99% level. The case for X-ray vs V band is 
identical. 



Energy bands Time lag [days] 

long term intensive 

X-ray vs B G-^IlIm S-^l^.s 

X-ray vs V 

B vs V 0.3^2.1 0At\-^2 



Table 1. Time delays in days between different pairs of energy 
bands, positive values indicate the first band leads. The lags were 
calculated using long-term, 4-day sampled light curves (middle 
column) and short-term, daily sampled light curves (right col- 
umn). 



curves. A first glance at Fig. |3] seems to indicate that the 
optical bands are leading the X-rays by ~ 10 days because 
of the sharp drop in flux in both bands around the middle of 
the light curve. The CCF analysis, however, reveals that the 
correlation at that lag only reaches a value of ~ 0.5 and that 
it is not significant. The only significant peak in the CCF 
is at positive lag values, with optical bands lagging, which 
make the overall intensive light curves match better. We 
note that the significant peak in the CCF at very small lag 
comes from the low amplitude X-ray/optical fluctuations, as 
shown in Fig. [S] 




V band lag [days] 

Figure 7. DCF derived using only the data from the period of 
intensive monitoring, shown in Fig. |3]in the X-ray and V bands. 
We used the daily-sampled V band data together with the 3- 
times daily sampled X-ray light curved, binned to 1 day bins, 
to estimate the central peak of the DCF more accurately. The 
central DCF peak is displaced to positive lags indicating that the 
short term V band fluctuations lag the X-rays by approximately 
5 days. The solid lines represent the 99% level for the correlation, 
estimated using simulated X-ray light curves of the appropriate 
length and sampling. The case for X-ray vs B band is identical. 



cance level curves, in this case we show the X-ray vs V band 
but note that the X-ray vs B band DCF is virtually identical. 
We note that the X-rays retain much more short time-scale 
variability power than the optical bands even after binning 
the X-ray data in 1-day bins. 

Although the light curves are significantly correlated, 
the CCF peak only reaches values ^^0.7, indicating that part 
of the variability in X-ray and optical bands is incoherent, 
which makes the delay harder to see directly from the light 



4.1 Time lags 

Time delays between the different light cures were esti- 
mated from the centroid of the central peak in the respec- 
tive correlation functions. The centroid was calculated as 
the weighted average of DCF values above 60% of the DCF 
maximum. The errors were calc ulated using the bootstrap 
method of lPeterson et all l|l998l ) by selecting randomly 67% 
of the points in each light curve and recording the lag cen- 
troid of the resulting DCF. We repeated this procedure 2000 
times for each pair of light curves to measure the median and 
spread of centroid values. The errors quoted correspond to 
the 67% limits of the distribution of centroid values on either 
side of the median. 

The time lags from the long term light curves are 
r = 6.6lg Q days between X-ray and B, r = S.Slg g between 
X-ray and V and r — 0.3^2 4 between B and V. The short, 
intensively sampled data produced lags of r = 5.7l2'g be- 
tween X-rays and B, r = G.Slg ]^ between X-rays and V, and 
T = 0.4lj^'2 between B and V. All the lag measurements are 
summarised in Table im positive lag values indicate higher 
energy band leading. In all cases the lags from the long and 
short term light curves are consistent. Both, V and B bands 
lag the X-rays significantly, from the intensively sampled 
light curves we measure these lags to be at least 3 and not 
more than 9 days long. The V and B band fluctuations are 
simultaneous within our time-resolution, with a delay of V 
behind B between -0.8 and 1.5 days. 

Figure |8] shows the B and X-ray light curves normalised 
to their mean fluxes, where the B data has been shifted back 
by 6 days, according to the time lag measured. It is clear that 
in this way the peaks in both light curves match very well, 
especially at low optical fluxes. The long term trend has 
slightly higher amplitude in the B band while the amplitude 
of the rapid fluctuations is much larger in the X-rays, as 
expected from the analysis of the PDS. 
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Figure 8. X-ray (black lines) and B band (blue dots joined by 
solid lines) light curves normalised to their mean fluxes. In the 
bottom panel the B light curve has been shifted back in time by 6 
days, i.e. by the time delay indicated in the DCF. The long term 
trends and short term (~ 10s of days) peaks are well aligned but 
the variability amplitudes are different. 



5 DISCUSSION 

The measured lag of optical behind X-ray fluctuations sug- 
gests a reprocessing scenario where at least part of the op- 
tical variability is produced by the variable X-ray heating. 
In this case, the delay, interpreted as a hght travel time, 
gives the distance between the X-ray corona and the loca- 
tion of the re processor. The black hole mass of NGC 3783 
is 3 X lO^M© (jPeterson et al.l[20o3 ) producing a hght cross- 
ing time for 1 gravitational radius of IRg/c = 150 s. The 
time delay of optical behind X-ray fluctuations of ~ 6 ± 3 
days therefore corresponds to a hght travel time through 
a distance of ~ 3500 ± 1750i?g. We note however that the 
region on the disc where the amount of reprocessed flux 
peaks will not necessarily correspond to the region of peak 
intrinsic emission due to internal heating, since X-ray heat- 
ing may dominate at larger radii as will be discussed below 
in Sec. 



5.1 Transfer functions in the reprocessing scenario 

Not all the optical variability can arise from reprocessed X- 
rays, however, as is evident from the power spectrum analy- 
sis. The PDS of the B band data shows that the variability 
power continues to increase towards lower frequencies, at 
least down to frequencies 4 x 10~* Hz, or correspondingly 
time-scales of ~ 300 days, while the X-ray variability power 
breaks and flattens for time-scales below 2 days. If all the 
optical variability was produced by reprocessing, the X-ray 
fluctuations on time-scales shorter than 100 days would have 
to be smoothed-out by differential light travel time to oppo- 
site sides of the reprocessor, requiring it to have a minimum 
radius of 50 hght days. In this case, the average delay be- 
tween optical and X-rays would be of the same order, much 
longer than the measured ~ 6 day lag. 

To exemplify this argument we constructed transfer 
functions for different reprocessing geometries. If the op- 



tical variability arises only from X-ray reprocessing then 
the transfer function has to produce both the smoothing 
of fluctuations down to time-scales of ~ 300 days and a 
delay of ~ 6 days between primary and reprocessed flux. 
In this case, the variable part of the optical hght curve 
would be o{t) = x{t) ® i{t), where x{t) is the X-ray hght 
curve, i{t) is the transfer function and ® denotes convolu- 
tion. The relation of their Fourier transform would therefore 
be 0(/) = X /(/), so that 



\o{f)\ 



(2) 



where |0(/)P and correspond to the PDS of the 

optical and X-ray light curves, respectively, as calculated in 
Sec. |3l The plots in Figure [9] compare this PDS ratio (cal- 
culated from the flts to the PDS using Eq. 1 and shown 
here in green triple-dot-dashed lines) to the func- 
tions predicted for each case. We note that the PDS ratio 
was calculated for optical and X-ray hght curves observed 
simultaneously, so the variations in the PDS shape expected 
due to the red noise nature of these hght curves is not rele- 
vant. 

The reprocessor geometries considered were a spherical 
shell distribution of small clouds, a flat optically thick disc 
and an optically thick disc which flares out at large radii. 
These were chosen to mimic structures that probably exist 
in the vicinity of the black hole such as an optically thick 
accretion disc and the broad line region (BLR) clouds. In all 
cases, model parameters were adjusted to produce the ob- 
served average delay of 6 days between primary and repro- 
cessed emission, calculated as the expectation value of the 
transfer function. The spherical shell produces the strongest 
smoothing of long-term trends for a fixed time lag, compared 
to the other geometries. This is because there is no delay be- 
tween the intrinsic emission and the response from the side 
of the sphere closest to the observer. In disc geometries (ex- 
cept for an edge-on view) the intrinsic flux must flrst travel 
to the nearest side of the reprocessor, which is located out of 
the line of sight, delaying the start of the response. There- 
fore, the same 'smoothing power' can be associated with a 
longer delay, e.g. if the disc is viewed face-on and the repro- 
cessing region is concentrated in a narrow ring then there 
will be no smoothing of the response but the delay can be ar- 
bitrarily long, for an arbitrarily large ring radius. Evidently, 
the spherical shell is the best bet to smooth out the long- 
term variability while producing a short lag, as required by 
the data, but we explore other geometries too for further 
analysis. We detail each model in the following subsections. 

5.1.1 Spherical shell 

The flrst geometry investigated comprised a spherically sym- 
metric distribution of small clouds that capture X-rays from 
the central source and re-emit a fraction of this flux in the 
optical bands. Given the symmetry of this setup the only pa- 
rameters are the radius of the shell R and the polar angle 9 of 
each position on the shell to the line of sight of the observer. 
The delay of the response with respect to the primary emis- 
sion is T{d) = R/c (1 — cos6'). Assuming that the clouds ra- 
diate isotropically and that they do not obscure one-another, 
the re-emitted flux seen by the observer at a given delay r is 
simply proportional to the area on the shell that corresponds 
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to that delay, i.e. dA(r)/dr = 2ttR^ smOde / dr = 2%R^, for 
r between and 2R/c. The that corresponds to this 

setup is shown by the solid line in the top panel of Fig. |9l 
where 7? was fitted to reproduce the observed lag. The 'rip- 
ples' in this profile are produced by the sharp drop of the 
transfer function at r = 2R/c. 

To smooth out this drop, the reprocessing region can 
span a range of radii rather than being a single thin shell. In 
this case the total transfer function will be the sum of the 
contributions from all sub-shells, weighted by the fraction 
of the sky covered by the clouds on each sub-shell C{R), 
i{t) = C{Rj)A{Rj). We used this setup to construct the 
rest of the functions shown in the top panel of Fig. O taking 
the thickness of each sub-shell as proportional to its radius. 

The functions shown in the top panel of Fig. [5] corre- 
spond to the thickest possible shells for each covering frac- 
tion distribution, C{R), using an inner radius of 10 Rg and 
outer radius fit to match the observed delay. The |/(/)|'^ 
functions changed smoothly from the thin shell profile to the 
corresponding thickest shell profile for intermediate configu- 
rations. The dashed line corresponds to a C{R) that peaked 
at the middle of the shell and dropped smoothly to zero at 
the inner and outer boundaries. The dot-dashed line corre- 
sponds to a constant C{R), where each sub-shell contributes 
an equal amount of reprocessed flux. The transfer function 
in this setup strongly peaks at small lags, where all sub- 
shells contribute and drops off gently towards larger lags. 
For very thick shells, as in the case shown in the figure, this 
sharp peak in the transfer function transfers some variability 
power on short time-scales, rising the high frequency tail of 
|/(/)|^ This effect is more extreme in the last case, shown 
by the dotted line, where the covering fraction drops with 
radius as 1/r so that thick shells make the transfer function 
even more peaked at small lags and most of the variability 
power is therefore transferred. We note that this configu- 
ration over predicts the short timescale variability without 
reducing the long term variability enough. 



6 I ' — ■ — ' — I 
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Figure 9. The ratio between the PDS model fits to the X-ray and 
B bands from Fig. 4 is shown by the triple-dot-dash lines in all 
panels. The other curves represent the squared Fourier transform 
of the transfer functions calculated for different geometries of a 
reprocessor, as labelled in each panel. For each geometry we cal- 
culated transfer functions for different sets of model parameters, 
with the constraint that the delay between primary and repro- 
cessed light curves should match the observed 6 days (see text 
for model and parameter descriptions). If the optical variability 
was entirely produced by reprocessing of the observed X-ray light 
curve then the PDS ratio would match the Fourier Transform of 
the transfer function of the correct reprocessor, i.e. Eq. 2 would 
hold. We see that the 6-day delay constraint strongly limits the 
smoothing power of the transfer function below ~ 2 X 10~^ Hz, so 
none of the reprocessors can smooth out optical variations enough 
while producing the observed delay. 



5.1.2 Truncated flat disc 

The second geometrical setup corresponded to the re- 
sponse of an infinite plane illuminated by an X-ray source 
at a height of lORg follow ing the prescription given in 
iKazanas fc NavakshinI l|200lf ). We varied the inner trunca- 
tion radius of the disc in order to fit the lag. The truncation 
is necessary because otherwise the reprocessed fiux peaks 
at a small radius, making the light travel time to this re- 
gion, and hence the delay, too short. Truncating the disc 
further out forces the reprocessing region to appear at a 
larger distance from the X-ray source and determines the 
delay . Intrinsic disc flux expecte d for a standard accretion 
disc (|Shakura fc SvunvaevI [l973l ) was added to the repro- 
cessed fiux, for each radius and time delay, to calculate the 
B band response. The curves shown in the middle plot of 
Fig. |9] correspond to a face-on view (least smoothing at long 
time-scales), an intermediate viewing angle of 45° and an 
edge-on view (most smoothing on long time-scales). The in- 
ner truncation radius needed to produce a 6 day lag in each 
case was 1950, 1550 and 1300 Rg, respectively. 
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5.1.3 Flared disc 

The third geometry corresponds to a flared disc, illuminated 
by a central X-ray source. The region of the disc where the 
scale height increases presents a large area facing the illumi- 
nating source and dominates the reprocessed emission. This 
scenario can therefore produce a large lag without requiring 
a large inner truncation radius of the disc. 

For the transfer function calculation, we assumed that 
the disc surface bends upwards sharply at a radius Rf and 
the normal to the flared surface makes an inclination angle 
a with respect to the plane of the flat section of the disc and 
has length 5*. We neglect differences in emissivity and time 
delays from different vertical positions on the flare surface. 
The effective area of this region that receives emission from 
a central primary source is dA = {Rf x S + sin a) cos ad6. 
If 9 is the azimuthal angle on the disc from the projection 
of the line of sight, the delay from each position is t{0) = 
Rf/c (1 — cos S sin i) and the area seen by the observer 
is dA'(T)/dT — {cosisina — sin i cos a cos 9)d A/ d6 dO/dr, 
where i is the viewing angle with respect to the axis of sym- 
metry of the disc. When i > a the near side of the disc is 
seen from behind, which makes dA' negative, as the disc is 
assumed to be optically thick, this area does not contribute 
to the transfer function. 

The curves shown in the bottom panel of Fig. [5] cor- 
respond to different viewing angles i and opening angles a. 
The strongest smoothing shown is produced for a = 70° and 
i = 45°, requiring a value of Rf = 3150-Rg. The interme- 
diate curves show cases where the opening and inclination 
angles are similar so the near side of the flared disc is viewed 
almost edge-on. In these cases the reprocessed flux concen- 
trates only at the far side of the disc, which is viewed more 
face-on. For this reason the response is retarded so the lag 
can increase while producing less long-term smoothing. The 
curves shown correspond to a = i = 75° and a = i = 45° 
with flare radii of Rf = 2350-Rg and Rf = 2700-Rg, respec- 
tively. The least amount of long-term smoothing is produced 
for face-on viewing angles, i = in this case shown with an 
opening angle of a = 45° and a radius Rf = 34OO-R9. 

5.2 Intrinsic plus reprocessed optical variability 

We note that in all the transfer functions described above, 
the required lag of 6 days strongly limits the smoothing 
power below ~ 2 x 10"'' Hz so the observed ratio of X- 
ray versus B band PDS cannot be reconstructed by a single 
transfer function. This adds to the fact that the B band 
has slightly larger long-term variations in normalised flux 
than the X-rays, which is hard to reconcile with a picture 
where these large amplitude optical trends are produced by 
reprocessing. We therefore conclude that another source of 
long term optical variability must exist and it is probably 
produced by the intrinsic optical emi tting region. 

As shown in lArevalo et al.l l|2008l ) rapid optical variabil- 
ity, imprinted by reprocessing on a long-term intrinsic fluc- 
tuation, can shift the CCF peak to the delay produced by 
reprocessing, regardless of the delay between the long term 
fluctuations. This is true even if the long term, intrinsic fluc- 
tuations in the optical bands have much larger amplitude 
than those produced by reprocessing. It is therefore possible 
that a small amount of the optical flux is produced by repro- 



cessing of X-rays at a large distance from the corona, while 
the rest of the optical emission and variability are produced 
intrinsically by the accretion disc at a different radius. In 
this case, the measured lag is the inprint of X-ray reprocess- 
ing on optically thick material at some distance from the 
X-ray source, which produces the rapid optical variability 
and 6-day delay. We now estimate the possible location of 
the optical intrinsic and reprocessed emitting regions. 

A radiatively efficient, optically thick disc radiates ap- 
proximately as a black body with a radially-dependent char- 
acteristic temperature and a total energy output defined by 
the gravitational energy loss. Assuming such a disc emits the 
optical flux observed in NGC 3783 we can calculate the loca- 
tion of the optical emitting regions. We use d the presc r iption 
of flux as a function of radius given in Treves et all (Il988l ) 
to calculate the local temperature and optical band flux for 
a mass of 3 x IO^Mq. We fitted the disc accretion rate in 
this formula to reproduce the observed optical flux obtain- 
ing a value of m = 0.01 of the Eddington value. We note 
that this accretion rate is lower than that estimated from 
the bolometric luminosity because it will only produce the 
thermal optical/IR emission from the accretion disc. It does 
not include the fraction of power directed into the corona 
(which results in X-ray emission) nor the additional IR flux 
directed into the line of sight by reprocessing structures such 
as the torus. The resulting 50% (90%) of the B band intrin- 
sic disc emission is contained within 73 (208) Rg. Given that 
the long term variability is probably intrinsic to this opti- 
cal emitting region, we estimate the relevant characteristic 
time-scales. At a radius of 100-Rg, which contains 60% of the 
B and 50% of the V flux, the dynamical time-scale (time it 
takes orbiting material to travel 1 rad of a Keplerian orbit) 
is 6.7 hours. Optical variability is observed on time-scales 
of 300 days, i.e. ~ 1000 times longer than the dynamical 
timescale at this radius, so accretion rate fluctuations on 
the local viscous time-scale are not an unreasonable source 
of this variabi lity. The viscous time-scale of a standard ct- 
disc (Shakura fcSvunvaevI 19731 ) is Tviac = Tdyn/{a{H / Rf), 
so in this case a{H/R)^ — 10"'^. Assuming a typical a — 0.1 
this implies a scale height-to-radius ratio H/R = 0.1, a rel- 
atively thick accretion disc. Alternatively, this variability 
time-scale could correspond t o the therm al time-scal e of the 
thin d isc, as discussed in e.g. iLiu et aL 

I ([2008, 1: Kellv et al.l 

(|2009l ). 

The location of the reprocessed optical emission depends 
on the location of the X-ray source and the geometry of the 
reprocessor. We constructed emissivity proflles for the B and 
V bands assuming a simple geometry of an accretion disc 
with constant H/R, and an X-ray source on the axis of sym- 
metry over the disc. The intrinsic disc accretion rate used 
was rh — 1% of the Eddington limit, a BH mass of 3 x 10^ Mq 
and a value of H/R = 0.1. The X-ray flux was equal to flve 
times the average 2-10 keV flux measured, to account for 
X-ray flux outside this band (see iDe Rosa et al. I l2002l , for 
their X-ray broadband analysis using BeppoSAX data) and 
assumed that 60% of the impinging X-ray flux were ther - 
malised by the disc (see discussion in lArevalo et al.l |2008| ) . 
Figure [10] shows the cumulative optical flux as a function 
of radius for the intrinsic, reprocessed and total emission 
for an X-ray source height of 50 Rg. The reprocessed curve 
was calculated as the difference between the intrinsic and 
total cumulative fluxes. Evidently, the intrinsic and repro- 
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Figure 10. Cumulative flux as a function of radius R in the B 
(solid lines) and V (dashed lines) bands. The intrinsic flux is emit- 
ted thermally by an accretion disc while an X-ray source added 
on the axis of symmetry enhances the disc emission by thermal 
reprocessing. The curves marked 'reprocessed' is the difference be- 
tween the intrinsic and total curves. The reprocessed flux peaks 
at a larger radius than the intrinsic emission and it is this larger 
value that must be used to interpret the delay between X-ray and 
optical fluctuations. 

cessed emissivities peak at different radii with the repro- 
cessed emission being produced further out. As this is the 
emission modulated by the variable X-ray flux, we should 
consider this outer radius to interpret the delay between X- 
ray and optical fluctuations produced by reprocessing. 

For a face-on view, the parameters used in Fig. [10] pro- 
duce a delay of ~ 10 hours for the X-rays to modulate up 
to 50% of the reprocessed flux in both B and V bands, too 
short to explain the observations. Using a thinner disc, down 
to H/R = 0.001 does not have a significant effect on this de- 
lay, nor does reducing the impinging X-ray flux down to the 
observed 2-10 keV value (instead of flve times this value as 
was used for the figure). Varying the position of the X-ray 
source does have a strong effect on the time delay, partly be- 
cause increasing the height moves the main reprocessing re- 
gion further out in the disc but mostly because of the longer 
light travel time down to the disc. The measured delay of 
a few days would require the X-ray source to be located at 
~ lOOO-Rj over the disc. 

A different geometry of the reprocessor, however, might 
also produce the correct time delay. For example a lower X- 
ray source and a flared disc would present a large area of the 
disc to the X-ray flux only at large radii, pushing the loca- 
tion of reprocessed emission further out. Alternatively, the 
reprocessed X-rays might come from a structure similar to 
the broad line region (BLR). Reverberation mapping cam- 
paigns of this source show that the broad H/3 line lags the 
UV contin uum variations by 1- 10 days, depending on the 
line width l|Peterson et al.ir2004l ). Therefore, the light travel 
time from the centre of the accretion disc (probable source of 
optical/UV continuum) to the BLR is similar to the distance 
between the X-ray corona and the reprocessing material. In 
this scenario, the broad line variability tracks the large am- 
plitude optical continuum fluctuations, which modulate the 
intrinsic optical flux originating close (--^ 100-Rg) to the black 
hole. At the same time, the lines can be modulated by the 



rapid optical fluctuations, produced by X-ray reprocessing 
at the same radius of the BLR, which already contains a 
delay to the intrinsic optical and X-ray emission. 

Recalling the three geometries of the reprocessor dis- 
cussed in Sec. 15.11 we can now judge their ability to pro- 
duce the rapid, small-amplitude optical fluctuations only. 
Although it is not straightforward to disentangle long and 
short term fluctuations, we can use the low flux section in 
the first half of the campaign to estimate the amplitude and 
timescales of what are probably signs of reprocessing. The 
two big fiares on this segment, which show a clear lag to 
the X-rays, have timescales of 10-20 days and amplitudes of 
a few mjy in the B band. Also, their relative amplitude is 
about half or less of the corresponding X-ray flares. The re- 
processing mechanism, therefore, should produce a few mJy 
of optical flux and not smooth out completely the variations 
on time-scales of tens of days. 

The truncated flat disc reprocessor seems unlikely, on 
energetics grounds because the X-ray source sees the disc 
almost edge-on, given the large truncation radius needed for 
the lag, for any reasonable source height. This small solid an- 
gle subtended by the disc produces too little reprocessed flux 
to make an observable signature in the optical light curves, 
by many orders of magnitude. This problem is circumvented 
when using the flared disc because the inclined or vertical 
portion of the disc can face the X-ray source directly, reaches 
a higher temperature and produces more optical flux. 

The BLR clouds can have large enough colu mn den- 
sities to be optically thick to X-rays (~ 10^'*cm~^. iNetzeil 
l2008i ). Reprocessing in the form of optical emission lines is 
possible and the effect of their ffuctuations on the B and V 
band ffuxes can be estimated from spectral monitoring. The 
re verberation mappin g campaign on NGC 3783 presented 
bv lStirpe et al. shows that the lines in this region of 

the spectrum contribute around 10% of the B band flux and 
vary by a maximum ratio of 30% in flux. Therefore, their 
contribution to the B and V band variability is too small to 
produce the observed rapid fluctuations of about 20%. The 
BLR, however, can also produce diffuse continuum emis- 
sion in addition to the lines, as shown by iKorista fc GoadI 
(|200ll ). These authors calculate the expected BLR contribu- 
tion to the optical/UV continuum in detail for the case of 
NGC 5548. They show that 10-20% of the observed B and 
V band continua can arise from BLR emission and that the 
majority of this flux is reprocessed rather than reflected. If a 
similar amount of continuum emission is contributed by the 
BLR in NGC 3783 and the X-rays (and possibly UV) inci- 
dent flux is modulating this emission, then the BLR clouds 
could correspond to the reprocessor that produces the rapid 
optical fluctuations and the 6-day delay. 

We calculated 'reprocessed' B band light curves by con- 
volving the observed X-rays with the transfer function of the 
flared disc and the spherical shell of Sec. 15.11 as these struc- 
tures are the most promising candidates for the reprocessor. 
We used the intensively sampled X-ray data to produce the 
simulated B band light curves and compare the resulting 
optical vs X-ray DCF with the DCF of the real data, to 
select the best flt in the r = ±30 day range. In the case 
of the flared disc we calculated B band emission assuming 
that the flared portion of the disc emits as a black body 
at the temperature corresponding to the reprocessed X-ray 
flux (as the radius is large, intrinsic heating at in this re- 
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Figure 11. DCF between observed B and X-ray bands in markers 
with error bars. The dashed line represents the DCF between the 
observed X-ray band and a simulated B band light curve produced 
by reprocessing the X-ray flux on a flared disc, the solid line 
corresponds to reprocessing on a spherical shell of clouds. The 
disc and sphere parameters were flt to match the observed DCF. 
Both reprocessors can produce the correct delay. The correlation 
peaks at a higher value for the simulated light curves, probably 
due to the additional sources of variability in the real optical light 
curve. 

gion is negligible) . As an additional constraint in the fitting 
routine we required an average reprocessed flux of at least 1 
mjy at the observer location, to reproduce the small ampli- 
tude optical fluctuations, which contain the 6-day delay. The 
height of the X-ray source was flxed at 10-Rg and the 2-10 
keV flux was amplifled by a factor of 5 as above, to include 
X-ray flux outside of this band. This setup could produce 
the required amount of B band flux and the correct delay 
for a range of parameter values. The best-fitting values were 
a = 60°, j = 45°, i?/ = 2200i?g and a length of the flared 
region of S = 600-Rg . In the case of the BLR reprocessor we 
only calculated the shape of the reprocessed light curve and 
relied on the numbers given by iKorista fc GoadI (|200lf l to 
produce sufficient reprocessed optical flux. The best-fitting 
geometrical parameters for this case where Rin — 1000i?g 
and a thickness of four times this radius. The fltted portion 
of the intensive DCF for both models is shown in Fig. 1111 
The model parameters are not well constrained in either case 
and the quality of the flts are similar so this analysis only 
shows that both structures are feasible reprocessors. 

The simulated reprocessed light curves lag the X-rays 
by approximately the same time as the real optical data, but 
their DCFs peak at a higher value than the DCF of the real 
data. We note that the simulated B band light curves are 
produced only through X-ray reprocessing, so a good corre- 
lation is expected. On the other hand, the real optical light 
curve probably contains other sources variability, producing 
the long term trends as argued above, which will lower the 
correlation peak. 



6 CONCLUSIONS 

We have monitored the Seyfert galaxy NGC 3783 over two 
years, simultaneously in the X-ray 2-10 keV and optical B 



and V bands. The flux in all bands is highly variable on all 
the time-scales probed. Our main conclusions drawn from 
the variability properties are detailed below. 

(i) We recalculated the power spectrum of NGC 3783 in 
the X-ray band. Our new, intensively sampled light curve 
was programmed to cover the time-scales where a possible 
QPO had been detected with earlier data ([Summons et al.l 
[2007). The high quality new data, however, shows that the 
entire PDS is consistent with a single-bend power law and 
the possible QPO is not confirmed. The PDS parameters we 
foun d are consistent with t he single-bend power law model 
fit of[s ummons et all (|2007l ). with values of = 0.8, an = 
2.2, = 5.8 X 10"® Hz. 

(ii) B and V bands vary simultaneously, with a delay of 
no more than 1.5 days. The light curves are almost identical 
but the fractional amplitude of variations is slightly larger in 
the B ba nd. These results are consistent with those obtained 
by Reichert et all (|l994, ) who monitored NGC 3783 in the 
UV and optical bands, finding the largest variability at the 
shortest wavelengths. 

(iii) The fluctuations in optical and X-ray bands are sig- 
nificantly correlated. The optical bands show a delay with 
respect to the X-rays of ~ 6 ± 3 days. Similar values are 
obtained for both B and V bands and for long-term and 
intensively sampled light curves. 

(iv) The PDS of the optical bands is consistent with a 
single power law of slope a = 1.6, where the variability 
power increases continuously toward the longest time-scale 
measured, ~ 300 days. Comparison with the PDS of the 
X-ray band shows that if all optical variability arises from 
reprocessing of X-rays, the reprocessor would have to be 
~ 50 light days across, to smooth out fluctuations on shorter 
time-scales. This size is inconsistent with the small value of 
the delay found between optical and X-ray fluctuations (6 
days) so another source of optical variability must exist. 

(v) We interpret the optical variability as arising from 
both reprocessing and intrinsic accretion rate fluctuations. 
Accretion rate variations on time-scales of hundreds of days 
can cause the large amplitude optical fluctuations and prop- 
agate to the X-ray emitting corona producing the long term 
correlation. Reprocessing of X-rays can produce the short- 
term, small amplitude optical variations explaining the mea- 
sured positive delay of optical with respect to X-rays. The 
most likely reprocessor geometries we investigated were an 
optically thick accretion disc which flares out at large radii 
and a spherical shell distribution of clouds. In the disc case 
the scale height increases at a radius similar to the loca- 
tion of the BLR, while in the second case the reprocessor 
could correspond to the BLR itself if the clouds can produce 
sufficient optical contin ua from reprocess e d X-r ays and UV 
photons, as discussed in lKorista fc GoadI l|200lh . 
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